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Abstract 


We  develop  two  new  models  for  gated,  pulsed  radio  frequency  interference  (RFI),  one  for  cycle 
slip  performance,  and  one  for  bit  error  performance.  We  compare  the  results  of  these  two  models  to 
an  existing,  approximate  model  that  is  currently  present  in  ITU-R  Report  M.2220,  herein  referred 
to  as  the  dynamic  duty  cycle  factor  (DDCF)  model.  After  deriving  our  new  models,  we  show 
through  case  studies  that  the  new  time- varying  models  often  predict  much  worse  performance  than 
the  corresponding  DDCF  counterparts.  We  therefore  conclude  that  further  investigation  (either  by 
simulation  or  laboratory  measurement)  is  warranted  to  verify  the  new  models  presented  herein. 

Keywords:  pulsed  RFI,  cycle  slips,  bit  errors,  duty  cycle. 
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Figure  1:  Graphical  depiction  of  pulsed  RFI  model.  The  red  box  depicts  a  window  of  time  over 
which  C/N0  is  averaged  to  produce  an  estimate  of  the  degraded  C/N0. 


1  Introduction 

Pulsed  Radio  Frequency  Interference  (pulsed  RFI)  is  a  phenomenon  that  can  degrade  the  perfor¬ 
mance  of  GPS  receivers  [1] — [9].  When  strong  pulses  that  lie  within  the  GPS  band  are  present  at 
the  input  of  a  GPS  receiver,  the  input  often  saturates,  which  has  the  effect  of  masking  any  useful 
data  during  the  time  periods  that  the  pulses  are  present.  The  effect  of  the  pulses  on  system  per¬ 
formance  is  a  function  of  the  receiver  type  (e.g.,  blanking  vs.  non-blanking,  fast  vs.  slow),  and, 
hence,  a  significant  amount  of  effort  has  been  put  into  modeling  the  effect  of  pulsed  RFI  on  various 
receiver  types. 

Recommendation  ITU-R  M.2030  [1]  provides  a  set  of  models  for  calculating  the  degradation  ef¬ 
fects  of  pulsed  RFI  by  computing  an  equivalent  continuous  wave  (CW)  noise  interference  that 
lowers  the  nominal  carrier-to-noise  density  C / No .  A  nominal  depiction  of  pulsed  RFI  interference 
that  is  considered  in  M.2030  is  shown  in  Fig.  1.  Here,  the  pulses  are  assumed  to  occur  with  a 
regular  pulse  repetition  frequency  (PRF).  The  post-correlator  effective  degradation  on  C/N0  for  a 
saturating  receiver  as  a  function  of  the  pulse  duty  cycle  PDC  is  [1]: 

AC /N0  =  201og(l  -  PDC)  (1) 

PDC  =  PW  x  PRF  (2) 

where  PW  represents  the  pulse  width  (sec)  and  PRF  is  the  pulse  repetition  frequency.  Eqn.  1 
represents  a  theoretical  model  of  how  to  estimate  C/N0  degradation,  while  the  red  box  in  Fig. 
1  shows  how  one  might  compute  this  degradation  via  measurement.  Instantaneous  C/No  values 
would  be  obtained  inside  the  window  and  the  results  would  be  averaged  to  obtain  a  crude  estimate 
of  C/No  when  the  pulsed  RFI  is  present.  If  the  pulsed  RFI  is  persistent,  it  is  easy  to  see  from  the 
figure  that  the  average  value  of  C/N0  is  roughly  independent  of  the  placement  of  the  window  and, 
therefore,  the  estimate  of  C/N0  is  an  approximately  time-invariant  quantity.  Therefore,  it  seems 
reasonable  to  use  this  model  of  degraded  C/No  as  an  indicator  for  the  four  basic  GPS  receiver 
functions  (namely  acquisition,  carrier  tracking,  code  tracking,  and  data  demodulation). 

NASA  has  proposed  a  plan  to  introduce  a  spaceborne  scatterometer  for  its  Solid  Moisture  Active 
Passive  (SMAP)  Program  that  will  be  operating  in  the  L2  frequency  band  [10].  GPS  receivers  op¬ 
erating  in  this  band  are  subject  to  corruption,  and  NASA  is  obliged  to  show  that  their  scatterometer 
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Figure  2:  Pulsed  RFI  waveform  with  gating.  The  pulses  are  present  for  robs  seconds  and  are  absent 
for  the  remaining  Ttc  -  robs  seconds. 


will  not  significantly  degrade  the  performance  of  a  GNSS  receiver  as  compared  to  ambient  opera¬ 
tion.  The  scatterometer  system  uses  a  mechanical  rotating  antenna,  which  can  illuminate  a  victim 
GPS  receiver  for  a  short  time  (jobs  =  300  msec)  during  each  rotation  period  (T TC  =  4.1  sec).  A 
graphical  depiction  of  this  new  type  of  pulsed  RFI  waveform,  where  the  pulsing  is  gated  to  model 
the  fact  that  the  high  frequency  pulses  will  only  be  present  at  the  receiver  for  300  msec,  is  shown 
in  Fig.  2.  There  is  a  single  sentence  in  ITU-R  Report  M.2220  which  states  that  the  corresponding 
degradation  in  C/N0  should  now  be  modeled  as 

AC/M)  =  201og(l  -  PDCum)  (3) 

PDCLIm  =  PDC^.  (4) 

The  degradation  model  of  Eqn.  3  and  4,  herein  referred  to  as  the  Dynamic  Duty  Cycle  Factor 
(DDCF)  model ,  states  that  one  can  view  the  gated  pulsed  RFI  waveform  as  having  a  composite 
duty  cycle  that  is  the  product  of  the  duty  cycle  of  the  pulsed  waveform  ( PDC )  and  the  duty  cycle 
of  the  of  the  gating  waveform  ( Tobs/TTc )•  While  such  a  model  is  clearly  convenient,  its  applicability 
to  modeling  the  real-world  effects  of  such  a  pulsing  waveform  is  questionable.  Effectively,  Eqn. 
3  states  that  the  effective  degradation  in  C/M0  can  be  computed  by  averaging  the  C/M0  values  in 
the  long  red  window  of  Fig.  2.  While  the  instantaneous  degradation  may  be  large  during  the  Tobs 
seconds  that  the  pulsing  is  present,  the  remaining  Ttc-  robs  seconds  of  quiescent  behavior  have  the 
potential  to  hide  the  negative  contributions  of  the  pulsing  waveform  when  the  long  term  average  is 
taken.  It  is  possible  — and  in  fact,  quite  likely — that  the  negative  effect  of  the  pulses  during  their 
"on"  time  can  impact  receiver  functionality  far  more  than  Eqn.  3  and  4  suggest. 

In  the  remainder  of  this  document,  we  develop  new  models  for  analyzing  the  impact  of  a  gated, 
pulsed  RFI  source  which  take  into  account  the  time-varying  nature  of  the  overall  pulsed  waveform. 
In  particular,  we  develop  one  model  which  allows  us  to  characterize  the  average  number  of  cycle 
slips  per  second  as  a  function  of  the  gated,  pulsed  RFI  parameters,  and  another  model  which  allows 
us  to  assess  the  probability  of  a  bit  error  in  data  demodulation  schemes.  We  compare  the  results 
of  these  models  to  the  corresponding  results  obtained  by  using  the  DDCF  degradation  model  of 
Eqn.  3  and  4,  and  we  effectively  show  that  the  time-varying  models  often  predict  much  worse 
performance  than  the  DDCF  counterparts.  We  present  our  results  for  a  variety  of  gating  duty 
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cycles,  pulsed  waveform  duty  cycles,  and  nominal  C/No  values  to  illustrate  the  large  range  of 
differences  between  the  two  models. 

2  Probabilistic  Model  for  Cycle  Slips 

Cycle  slips  are  an  important  measure  of  GPS  receiver  performance  [5].  They  occur  when  the 
phase  estimate  of  the  carrier  waveform  shifts  by  360°.  This  offset  is  important  as  it  degrades  the 
accuracy  of  a  PNT  solution.  Moreover,  when  the  incoming  signal  is  encoded  with  data,  cycle  slips 
(now  of  180°)  can  cause  data  bit  inversions,  thus  creating  bit  errors.  Here  we  present  a  model  for 
assessing  the  average  number  of  cycle  slips  per  second  (ACS)  as  a  means  of  quantifying  cycle  slip 
performance,  which,  as  we  shall  see,  is  largely  linked  to  C/No. 

When  C/No  is  constant,  the  number  of  cycle  slips  k  in  a  a  period  of  time  t  is  Poisson  distributed 
[4,  5].  Specifically,  if  we  denote  by  Pr{k(t )}  the  probability  of  k  cycle  slips  occurring  in  t  seconds, 
we  have 
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where  TCs  is  called  the  mean  time  to  cycle  slip ,  and  the  parameters  in  Eqn.  6  are  as  listed  in  the 
table  below: 


•  6e\  steady-state  jerk  angle. 

•  Bl :  one  sided  loop  bandwidth. 

•  C  third  order  implementation  loss. 

•  p:  period  of  equivalent  S-curve  of  tracking  loop. 

•  erf  RMS  phase  error  of  tracking  loop. 

•  /„():  77-th  order  modified  Bessel  function  of  the  first  kind. 


For  carrier  tracking,  the  parameter  <xf/)  above  is  given  by 
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Figure  3:  Graphical  depiction  of  periodically  varying  C/No  for  a  low  frequency  gating  waveform. 


In  Eqn.  7  and  8,  TL  refers  to  the  coherent  integration  time,  and  crR  and  crT  refer  to  the  receiver  and 
satellite  RMS  noise  levels,  respectively.  It  is  clear  from  Eqn.  8  that  the  RMS  phase  error  is  highly 
a  function  of  C/No.  In  particular,  as  C/N0  decreases,  <t(!)  increases,  and  Tcs  decreases,  meaning 
that  cycle  slips  become  more  frequent. 


By  virtue  of  the  Poisson  distribution  in  Eqn.  5,  the  expected  number  of  cycle  slips  in  a  given  time 
E{k(t)}  is 

E{k(t)}  =  (9) 

I  cs 


Observe  that  the  expected  number  of  cycle  slips  grows  linearly  with  time.  We  define  ACS,  the 
average  number  of  cycle  slips  per  second,  as  the  long-term  average  of  the  expected  number  of 
cycle  slips: 


ACS 


a  r  E{k(t)} 
—  lim - 


(10) 


When  C/No  is  constant,  Eqn.  10  tell  us  that  the  ACS  is  simply  the  reciprocal  of  the  mean-time-to- 
cycle  slip,  an  intuitive  result. 


2.1  ACS  for  Periodically  Varying  C/N{) 

We  model  the  effect  of  gating  on  the  high  frequency  pulsed  waveform  by  considering  C/No  that  is 
periodically  varying,  as  depicted  in  Fig.  3.  When  pulsed  RFI  is  not  present  (corresponding  to  the 
case  when  the  rotating  antenna  is  not  illuminating  the  GPS  receiver),  then  the  value  of  C/N0  is  some 
nominal  value  (C/M))nom.  When  the  pulsed  RFI  is  present,  the  effective  noise  level  increases,  and 
C/No  is  lowered  for  T \  seconds  to  a  value  we  denote  as  (C/No)rfi-  This  process  repeats  periodically 
with  period  T2.  In  terms  of  the  NASA  SMAP  example  discussed  in  the  introduction,  T i  =  and 
T2  =  Ttc  (refer  to  Eqn.  4). 

To  derive  an  expression  for  the  ACS  in  this  new  scenario,  we  first  need  to  model  the  probability 
that  k  cycle  slips  occur  in  a  given  time  t.  Because  C/No  now  varies  with  time,  the  simple  Poisson 
model  of  Eqn.  5  does  not  directly  apply.  To  begin,  let  us  first  find  the  probability  mass  function 
for  the  number  of  cycle  slips  k  that  occur  in  one  period  of  the  C/No  waveform  T2.  To  aid  us  in 
our  computation,  it  is  helpful  to  partition  the  time  length  T2  into  two  separate  regions:  a  region 
of  length  T\  where  C/N(]  =  (C//V0)rfi  and  a  region  of  length  T2-T\  where  C/No  -  (C/No) nom- 
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Observe  now  that  the  following  equality  holds 


k 

Pr {k{T2)}  =  Pr{/n(T1)}Pr{(&  -  m){T1  -  Ty)}  (1 1) 

772=0 


The  above  relationship  says  the  following:  suppose  that  a  total  of  cycle  slips  that  occurs  in  time 
T2  is  k,  and  that  of  these  k  cycle  slips,  in  of  them  occur  in  the  first  T \  seconds.  It  follows  then  that 
k  -  m  cycle  slips  occur  in  the  remaining  time  T2-T\.  The  probability  that  k  cycle  slips  occurs  in 
time  T2  is  therefore  the  sum  over  all  possible  combinations  of  having  in  cycle  slips  occur  in  time 
T i  and  k  -  m  cycle  slips  occur  in  T2.  Observe  that  because  C/N0  is  constant  over  each  interval: 


Pr  {m(Ty)} 

Pr  {(k-m)(T2-Ty)} 


(k  -  in) ! 


(12) 
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Substituting  Eqn.  12  and  13  into  Eqn.  11  and  evaluating  the  sum  yields 


Pr  {k(T2)} 
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Hk exp  (~n) 


k\ 


T 1  T2  -  Ty 

— -  +  — — - 

Tcs  RFI  Tcs  nom 


(14) 

(15) 


The  above  result  is  likely  familiar  to  those  who  know  queuing  theory,  as  it  is  simply  a  restatement 
that  independent  Poisson  processes  regenerate  under  addition.  Those  seeking  a  proof  of  this  state¬ 
ment  are  referred  the  Appendix.  While  the  probability  mass  function  (PMF)  for  the  number  of 
cycle  slips  is  not  Poisson  for  an  arbitrary  time  t,  it  behaves  as  a  Poisson  random  variable  when 
the  time  interval  under  consideration  is  the  period  length  T2.  Moreover,  it  can  be  shown  that  for 
integer  multiples  n  of  the  period  length  T2  (see  the  Appendix  for  a  proof): 


Pr  {k(nT2)} 
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Hk  exp(-/r) 
k\ 

nTy  +  n(T2  -Ty) 

TcS  RFI  TcS  nom 


(16) 

(17) 


A  direct  result  of  Eqn.  16  is  that  E{k(nT2)}  =  n  of  Eqn.  17.  We  can  use  this  result  to  compute  the 
ACS  for  the  time-varying  C/N()  waveform.  We  are  interested  in  computing 


ACS  =  lim 

7— >oo 


E{k(t)} 

t 


(18) 


To  compute  this  limit,  we  first  note  that  the  E{k{t)}  is  monotonically  increasing  with  t.  Indeed, 
because  the  arrival  times  of  individual  cycle  slips  are  modeled  as  independent,  then  whenever 


t2  >  ty, 


E{k(t2)}  =  E{k(ty)}  +  E{k(t2  -  ty)}  >  E{k{ty)}. 


(19) 
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In  particular,  whenever  nT2  <  t  <  (n+  1  )T2,  it  follows  that  E{k(nT2)}  <  E{k(t)}  <  E{k((n  +  1)T2)} 
and,  moreover,  that 


E{k(nT2)}  <  E{k(t)}  <  E{k((n  +  \)T2)} 


in  +  DTt 


nT2 


By  the  squeeze  theorem, 


E{k(nT2)}  E{k(t)}  E{k  ((n  +  l)T2)} 

lim  - - <  lim - <  lim  — 


( n  +  \)T2  t 

Utilizing  Eqn.  17,  we  find  that 

r  E{k(nT2)}  E{k((n  +  1)T2)} 

lim - =  lim - = 

n—>oo  (j'l  +  1  )T2  n — »oo  nT2 

which,  by  virtue  of  the  inequality  chain  in  Eqn.  21,  leads  to 
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The  result  of  Eqn.  23  has  a  somewhat  intuitive  appeal.  If  we  define  the  variables 

A  1 

\CK  —  _ 

u  nnm  —  — 


ACS  rf,  ^ 
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1 

CS  RFI 


(23) 


(24) 

(25) 


then  Eqn.  23  is  equivalent  to 


ACSeff  =  ^ACSRFI  + (l-^j 


ACS, 


(26) 


where  the  "eff"  subscript  is  used  to  denote  the  effective  ACS  of  the  time-varying  model.  Eqn. 
26  states  that  the  effective  ACS  when  C/Mi  varies  periodically  is  a  weighted  average  of  two  time- 
invariant  ACS’s,  one  corresponding  to  the  nominal  C/N0  and  one  corresponding  to  the  C/N()  during 
the  pulsed  interference. 


2.2  Procedure  for  Comparing  ACS:  Time-varying  Model  vs.  Av¬ 
erage  C/Nq  Model 

We  now  describe  the  method  for  computing  ACSeff  using  our  time-varying  model  and  compare  it 
to  the  ACS  ddcf?  the  ACS  that  is  predicted  by  the  DDCF  model  of  Eqn.  3  and  4.  After  describing 
the  general  methods  in  both  cases,  we  shall  present  results  for  a  variety  of  case  studies. 
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2.2.1  Computing  ACS  Using  Time-varying  Model  (Eqn.  26) 


As  mentioned  in  the  introduction,  ITU-R  Report  M.2030  provides  a  model  for  computing  the 
degradation  in  C/No  when  non-gated,  pulsed  RFI  is  present  (repeated  here  for  convenience): 


A  C/No  =  201og(l  -PDC) 


(27) 


where  PDC  is  the  duty  cycle  of  the  high  frequency,  pulsed  RFI  signal.  To  compute  ACS  using  the 
time-varying  model  of  Eqn.  26,  we  use  the  following  procedure: 

•  For  a  nominal  value  of  (C/M0)n om  (when  the  pulsing  is  turned  off  by  the  gating  waveform), 
compute  TCSnom  using  Eqn.  6. 

•  Use  the  computed  value  of  TCSnom  to  compute  ACSnom  via  Eqn.  10. 

•  For  a  given  value  of  PDC,  compute  (C/MOrh  =  (C/M0)no m  +  A C/N0,  where  AC /No  is  com¬ 
puted  via  Eqn.  27. 

•  Use  (C/A0)rfi  to  compute  ?CsRFI  via  Eqn.  6  and  then  ACS  RFi  via  Eqn.  10. 

•  For  given  values  of  T\  =  robs  and  T2  =  Ttc,  compute  ACAeff  via  Eqn.  26. 


2.2.2  Computing  ACS  Using  the  DDCF  Model  (Eqn.  3  and  4) 

Also  mentioned  in  the  introduction,  ITU-R  Report  M.2220  [2]  and  ITU-R  Recommendation  M.2030 
[1]  provide  a  model  for  computing  degradation  in  C/N0  when  gated,  pulsed  RFI  is  present  (repeated 
here  for  convenience): 


A  C/No  =  20  log(l  -PDCum) 


(28) 

(29) 


where  PDC  is  the  duty  cycle  of  the  high  frequency  pulsed  RFI  waveform,  T TC  is  the  period  of  the 
gating  waveform,  and  Tobs  is  the  amount  of  time  during  one  gating  period  for  which  the  pulsed  RFI 
is  present.  To  compute  ACS  using  the  DDCF  model  (denoted  ACS  ddcf),  we  use  the  following 
procedure: 

•  For  given  values  of  PDC,  robs  and  Ttc,  compute  (C/N0) ddcf  =  (C/N0)nom  -  AC/M),  where 
AC/M)  is  computed  via  Eqn.  28  and  29. 

•  Using  the  calculated  value  of  ( C/M0)ddcf ,  compute  TCsddcf  via  Eqn.  6. 

•  Using  the  calculated  value  of  7CSddcf,  compute  ACS  ddcf  via  Eqn.  10. 
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Parameter 

Value 

0e 

o 

o 

Bl 

10Hz 

s 

-l.OdB 

p 

2n,  Pilot  Tone 

7 r,  with  Data 

Tl 

20msec 

CTr 

2.9° 

(T  j 

2.2° 

Table  1:  Parameter  values  used  in  ACS  calculations  that  were  used  to  generate  Fig.  4  through  Fig. 
7. 

2.2.3  Results 


Figures  4  through  7  show  the  results  of  computing  the  ACS  via  the  two  different  models  for  a 
variety  of  scenarios.  In  each  plot,  the  horizontal  axis  represents  the  duty  cycle  of  the  pulsed  wave¬ 
form  from  1%  to  50%,  and  the  vertical  axis  shows  the  corresponding  ACS  (the  gating  duty  cycle  is 
held  constant  for  each  plot).  The  dashed  red  waveform  indicates  the  ACS  predicted  via  the  DDCF 
model,  ACS  ddcf,  whereas  the  blue  waveform  represents  the  ACS  predicted  by  the  time- varying 
model,  ACSeff.  The  nominal  value  of  C/No  is  shown  in  the  caption  for  each  plot,  and  the  values 
for  the  various  parameters  that  are  inputs  to  Eqn.  6  for  computing  Tcs  are  shown  in  Table  1. 

Two  different  signal  scenarios  are  depicted  in  the  figures.  Fig.  4  and  Fig.  5  show  ACS  results  for 
the  case  where  the  signal  being  received  is  a  pilot  tone  (no  data  modulation).  Such  a  scenario  is 
characteristic  of  L2C  receivers.  In  this  case,  a  cycle  slip  occurs  when  a  full  2/r  radians  of  phase 
error  has  accumulated,  which  accounts  for  the  choice  of  p  =  2/r  in  Table  1.  By  contrast,  Fig.  6 
and  Fig.  7  depict  ACS  results  when  the  received  signal  has  modulated  data.  In  this  case,  a  cycle 
slip  can  occur  only  after  n  radians  of  phase  error  (due  to  ambiguity  of  the  data  bit  multiplication  of 
±1). 

Some  comments  are  in  order.  First,  observe  that  whenever  the  gating  duty  cycle  is  100%,  the  two 
models  produce  the  exact  same  result.  Such  is  to  be  expected  since  the  original  (non-gated)  pulsed 
RFI  waveform  is  obtained  whenever  gating  is  100%.  Also,  observe  that  whenever  PDC  is  small, 
the  two  models  also  effectively  produce  the  same  results  (when  PDC  is  zero,  there  is  no  pulsed  RFI, 
and  the  two  models  should  again  produce  identical  results).  As  PDC  increases  from  zero,  the  two 
models  produce  strikingly  different  results,  with  the  time-varying  model  producing  significantly 
higher  estimates  AC5efj  than  that  of  ACS  than  the  DDCF  model  ACS  ddcf-  This  effect  is  most 
notable  for  small  gating  duty  cycles  (e.g.,  7%),  but  is  also  pronounced  for  moderate  gating  duty 
cycles  (e.g.,  25%)  as  well. 

Overall,  the  depicted  results  indicate  that  the  DDCF  model  may  be  overly  optimistic  in  predicting 
ACS.  While  some  simulated  and/or  measured  data  should  be  obtained  to  verify  the  time-varying 
model  developed  here,  at  the  very  least,  the  results  suggest  that  one  should  use  extreme  caution  in 
applying  the  DDCF  model  for  PDC  values  that  are  more  than  a  few  percent. 
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Figure  4:  ACS es  vs.  ACS Ddcf  for  a  pilot  tone,  (C/No)mm  =  30dB 


Av^age  cycle  slips/sec,  Nominal  C/NO  =  30dB,  Gating  Duty  Cycle  =  7%,  Pilot  Tone 


(a)  Pilot  Tone,  Gating  Duty  Cycle  =  7% 


Av^age  cycle  slips/sec,  Nominal  C/NO  =  30dB,  Gating  Duty  Cycle  =  25%,  Pilot  Tone 


(b)  Pilot  Tone,  Gating  Duty  Cycle  =  25% 


Average  cycle  slips/sec,  Nominal  C/NO  =  30dB,  Gating  Duty  Cycle  =  100%,  Pilot  Tone 


(c)  Pilot  Tone,  Gating  Duty  Cycle  =  100% 
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Figure  5:  /4CS'eir  vs.  ACS Ddcf  for  a  pilot  tone,  (C/No)mm  =  25dB 


(a)  Pilot  Tone,  Gating  Duty  Cycle  =  7% 


Average  cycle  slips/sec,  Nominal  C/NO  =  25dB,  Gating  Duty  Cycle  =  25%,  Pilot  Tone 


(b)  Pilot  Tone,  Gating  Duty  Cycle  =  25% 


Average  cycle  slips/sec,  Nominal  C/NO  =  25dB,  Gating  Duty  Cycle  =  100%,  Pilot  Tone 
10 
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10 
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(c)  Pilot  Tone,  Gating  Duty  Cycle  =  100% 
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Figure  6:  /\CS'eir  vs.  ACS  ddcf  for  a  data-modulated  signal,  (C/No) nom  =  30dB 


Average  cycle  slips/sec,  Nominal  C/NO  =  30dB,  Gating  Duty  Cycle  =  7%,  with  Data 


10 

10 

a 

S  10 

CO 

a. 

"55 

*  10 

o 

o 

$  10 
< 

10 

10 

10 


(a)  Data-modulated,  Gating  Duty  Cycle  =  7% 


Average  cycle  slips/sec,  Nominal  C/N0  =  30dB,  Gating  Duty  Cycle  =  25%,  with  Data 
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(b)  Data-modulated,  Gating  Duty  Cycle  =  25% 


(c)  Data-modulated,  Gating  Duty  Cycle  =  100% 
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Figure  7:  ACS efr  vs.  ACSddcf  for  a  data-modulated  signal,  (C/No) nom  =  25dB 


(a)  Data-modulated,  Gating  Duty  Cycle  =  7% 


(b)  Data-modulated,  Gating  Duty  Cycle  =  25% 


(c)  Data-modulated,  Gating  Duty  Cycle  =  100% 
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3  Bit  Error  Probability  for  LI  C/A  and  P(Y)  Codes 


We  now  turn  our  attention  to  modeling  the  effect  of  gated,  pulsed  RFI  on  data  bit  errors.  The 
cycle  slip  analysis  of  the  previous  section  provides  us  some  insight  into  how  time- varying  C/N0 
will  affect  the  probability  of  bit  errors  (cycle  slips  will  cause  data  bit  inversions  on  data  modulated 
channels).  Here  we  perform  a  direct  analysis  of  data  bit  error  probability  for  a  specific  scenario, 
namely  LI  C/A  (or  P(Y)  code)  which  uses  a  (32,26)  Hamming  code  to  decode  the  data.  In  this 
particular  case,  for  constant  C/No,  the  probability  of  a  bit  error  conditioned  on  knowledge  of  carrier 
phase  error  (pe  is  given  by  [5] 


-erfc 

2 


(C/Np)  cos2  </>t 
Rn 


(30) 


where  Rb  is  the  bit  transmission  rate.  The  carrier  phase  error  (be  is  itself  a  random  variable  with 
PDF  [5], 

^  cos{^)e— 6 e)^ 


exp 


P(</>e)  = 


2  nlr 1  1 


0  Z2 


\(pe\  <  n- 


(31) 


where,  again,  a1,  is  computed  via  Eqn.  7  and  8.  The  average  probability  of  bit  error  is  thus 


Pb  = 


PB\^eP{(f>e)d(f>e 


(32) 


When  C /Nq  is  piecewise  constant,  as  in  Fig.  3,  we  can  compute  the  average  probability  of  error 
using  Bayes’  Rule.  Suppose  we  denote  by  the  average  probability  of  a  bit  error  when  C/No  = 
(C/No)nom,  and  P^FI  as  the  average  probability  of  a  bit  error  when  C/N0  =  (C/MOrfi-  Then  the 
effective  average  probability  of  a  bit  error  /J/fl  is  given  by 


Pf'P{C/N0  =  (C/A0)rfi}  +  PnBomP{C/No  =  (C/jVo)nom} 


(33) 


3.1  Procedure  for  Comparing  Average  Bit  Error  Probability:  Time- 
varying  Model  vs.  Average  C/N0  model 

The  procedure  we  use  to  compare  the  average  bit  error  probability  obtained  via  the  time- varying 
model  to  that  obtained  by  the  DDCF  model  is  very  similar  to  the  procedure  for  comparing  ACS 
for  the  two  models.  We  describe  the  exact  procedures  here  and  then  present  results  for  a  variety  of 
case  studies. 

3.1.1  Computing  Average  Probability  of  Bit  Error  for  Time-varying  Model 
(Eqn.  34) 

•  For  a  nomminal  value  of  (C/N0)nom  (when  the  pulsing  is  turned  off  by  the  gating  waveform), 
compute  o-0  via  Eqn.  7  and  8. 
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Parameter 

Value 

0e 

0.35° 

Bl 

10Hz 

Tl 

10msec 

O-R 

2.9° 

cr  j 

2.2° 

Rb 

50  bits/sec 

Table  2:  Parameter  values  used  in  PB  calculations  that  were  used  to  generate  Fig.  8  and  Fig.  9. 

•  Substitute  the  computed  value  of  cr 0  into  Eqn.  30  and  use  this  to  compute  via  Eqn.  32. 

•  For  a  given  value  of  PDC,  compute  (C//V0)rfi  =  (C/N0)nom  +  AC/N0,  where  AC /No  is  com¬ 
puted  via  Eqn.  27. 

•  Use  (C/No)rfi  to  compute  the  corresponding  value  of  (T(b  via  Eqn.  7  and  8,  and  use  this  value 
to  compute  P^F1  via  Eqn.  32. 

•  For  given  values  of  =  robs  and  T2  =  Ttc,  compute  Pe^  via  Eqn.  34. 


3.1.2  Computing  Average  Bit  Error  Probability  Using  DDCF  Model  (Eqn.  3 
and  4) 

The  procedure  for  computing  average  bit  error  probability  using  the  DDCF  model,  denoted  by 
PfCF),  is  again  very  similar  to  the  procedure  for  computing  ACS  using  the  DDCF  model: 


•  For  given  values  of  PDC,  robs  and  Ttc,  compute  (C/M0) ddcf  =  (C/Mdnom  +  AC/M),  where 
AC/M)  is  computed  via  Eqn.  28  and  29. 

•  Using  the  calculated  value  of  ( C/M0)ddcf ,  compute  the  corresponding  value  of  <x(4  via  Eqn. 
7  and  8. 

•  Substitute  (Tlb  into  Eqn.  30,  and  use  this  expression  to  compute  PgDCF  via  Eqn.  34. 

3.2  Results 

Figures  8  and  9  depict  the  results  of  computing  the  average  bit  error  probability  for  two  different 
scenarios.  As  before,  the  horizontal  axis  in  each  plot  represents  the  duty  cycle  of  the  pulsed  wave¬ 
form  from  1%  to  50%,  and  the  vertical  axis  shows  the  corresponding  average  bit  error  probability. 
The  dashed  red  waveform  depicts  PgDCF  while  the  blue  waveform  depicts  P^  for  the  time- varying 
model.  The  nominal  values  of  C/N0  are  shown  in  the  caption  for  each  plot,  and  the  values  for 
various  parameters  that  inputs  to  the  relevant  equations  are  shown  in  Table  2. 

The  plots  in  Fig.  8  and  Fig.  9  depict  the  same  qualitative  behavior  as  their  ACS  counterparts 
Fig.  4  through  Fig.  7;  the  two  models  coincide  for  gating  duty  cycles  of  100%  and  for  pulse  duty 
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cycles  of  0%.  Furthermore,  for  low  gating  duty  cycles,  quickly  grows  to  be  much  larger 
than  PgDCF  as  PDC  increases.  Therefore,  the  same  general  conclusions  and  warnings  apply  to  the 
average  probability  of  bit  error  as  for  ACS:  the  DDCF  model  exhibits  overly  optimistic  behavior, 
so  extreme  caution  should  be  used  when  applying  it  for  situations  where  the  PDC  values  are  more 
than  a  few  percent. 

4  Summary 

Our  goal  in  this  document  has  been  to  evaluate  the  validity  of  the  DDCF  model  that  is  currently 
present  in  ITU-R  Recommendation  M.2030  [1].  Using  the  average  number  of  cycle  slips  per  sec¬ 
ond  and  the  average  probability  of  a  bit  error  as  high-level  performance  indicators,  we  developed 
models  which  explicitly  take  into  account  the  time- varying  nature  of  the  C/No  profile  and  com¬ 
pared  their  performance  to  the  approximate  DDCF  model.  Our  results  indicate  that  in  several 
realistic  scenarios,  the  newly  developed  time-varying  models  predict  much  worse  performance 
than  the  corresponding  DDCF  model  and,  therefore,  the  applicability  of  the  DDCF  model  should 
be  questioned.  Further  analysis,  either  via  simulation  or  laboratory  measurement,  is  warranted  to 
assess  the  validity  of  the  time-varying  models  presented  herein. 
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Figure  8:  Pf  vs.  P°DCF,  (C/N0) nom  =  30dB 


(a)  Gating  Duty  Cycle  =  7% 


(b)  Gating  Duty  Cycle  =  25% 


(c)  Gating  Duty  Cycle  =  100% 
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Figure  9:  Pf  vs.  P°DCF,  (C/N0) nom  =  25dB 


(a)  Gating  Duty  Cycle  =  7% 


(b)  Gating  Duty  Cycle  =  25% 


(c)  Gating  Duty  Cycle  =  100% 
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Appendix  A 


A.1  Proof  of  Eqn.  12  and  13 


To  simplify  notation,  we  make  the  following  definitions: 

T\ 


A 

a  — 


P  = 


TcStm 
TWi 
T cs 

nom 


With  the  above  substitutions,  the  sum  of  Eqn.  11  becomes 

a,ne-a  pk-me-P 
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(34) 

(35) 


(36) 

(37) 

(38) 

(39) 

(40) 


Observe  that  the  binomial  theorem  was  applied  in  Eqn.  38  to  obtain  Eqn.  39.  Noting  that  n  =  a+fi 
completes  the  proof. 


Implicit  in  the  above  derivation  is  the  assumption  that  we  consider  a  time  interval  of  length  T2 
which  starts  on  the  falling  edge  of  the  C/N0  waveform  of  Fig.  3  and  ends  at  the  next  falling  edge. 
This  assumption,  however,  is  not  critical  and  the  result  holds  for  an  arbitrary  interval  of  length  T2 
as  we  now  show.  Suppose  we  consider  an  interval  of  length  T2  which  starts  t  units  of  time  before  a 
falling  edge  of  the  C/No  waveform.  Then  the  probability  that  k  cycle  slips  occur  in  an  interval  that 
ends  T2  seconds  after  this  starting  time  can  be  expressed  as 


Pr  {k(T2)}  = 
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(41) 

(42) 

(43) 

(44) 
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Observe  now  that 


Pr  {k(T2)} 


k  (  JIL,  ane~a  y0"~'"V_'8  ' 
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(45) 

(46) 

(47) 


In  the  above,  we  have  applied  the  first  part  of  this  proof  twice  in  evaluating  the  sums.  Noting  that 
a  +/3  +  y  =  /j.,  we  see  that  Eqn.  47  is  equivalent  to  Eqn.  14.  A  similar  set  of  sums  can  be  written  for 
the  situation  that  the  initial  starting  point  for  the  length  T2  interval  lies  between  the  falling  edge  of 
the  C/Nq  waveform  and  the  rising  edge,  and  the  result  is  exactly  the  same.  Hence,  Eqn.  14  and  15 
represent  the  PMF  for  the  number  of  cycle  slips  that  occur  in  an  arbitrary  time  interval  of  length 
T2. 


A.2  Proof  of  Eqn.  16  and  17 

The  proof  of  this  statement  is  very  similar  to  the  last  proof  of  Eqn.  14  and  15  and,  hence,  many  of 
the  details  are  omitted  for  brevity.  Note  that  we  have  already  proved  the  base  case  n  =  1 .  Assuming 
the  statement  is  true  for  some  arbitrary  value  of  n,  then  the  probability  that  k  cycle  slips  occur  in 
time  ( n  +  1  )T2  can  be  expressed  as 

k 

Pr {k ((n  +  1)T2)}  =  Y,  prW  (nT2)\Pv\(k  -  m)  (T2)|  (48) 

m= 0 

Stated  in  words,  the  probability  that  k  cycle  slips  occur  in  time  ( n  +  1  )T2  is  the  probability  that  m 
cycle  slips  occur  over  the  first  nT2  seconds  followed  by  k  -  m  cycle  slips  occur  in  the  following  T2 
seconds,  summed  over  all  possible  values  of  m.  We  have  already  established  that  the  two  quantities 
in  the  sum  in  Eqn.  48  are  Poisson  random  variables.  Hence,  the  same  computation  may  be  carried 
out  as  in  the  Proof  of  Eqn.  14  and  15  to  show  the  general  result. 
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